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ABSTRACT 

Using results from high-resolution galaxy formation simulations in a standard ACDM cosmology 
and a fully conservative multi-resolution radiative transfer code around point sources, we compute 
the energy-dependent escape fraction / esc of ionizing photons from a large number of star forming 
regions in two galaxies at five different redshifts from z = 3.8 to 2.39. All escape fractions show a 
monotonic decline with time, from (at the Lyman-limit) ^6 — 10% at z = 3.6 to ~ 1 — 2% at z = 2.39, 
due to higher gas clumping at lower redshifts. It appears that increased feedback can lead to higher 
/ csc at z > 3.4 via evacuation of gas from the vicinity of star forming regions and to lower / csc at 
z ^ 2.39 through accumulati on of swept-up shel ls in denser environments. Our results agree well with 
the observational findings of llnoue et alJ l{2006[) on redshift evolution of f csc in the redshift interval 
2=2-3.6. 

Subject headings: galaxies: formation — intergalactic medium — HII regions — radiative transfer 



1. INTRODUCTION 

In the redshift interval 2 < z < 6 most ionizing photons 
in the Universe are thought to originate in PopII stars in 
norm al (V c > 30 — 50kms _1 ) galaxies ijNagamine et al.l 
120061) . For reionization calculations at z > 6 photons 
from lower mass halos should be taken into account 
ijlliev et al.ll2006l) , whereas at z ~ 2 active galactic nuclei 
start to play a dominant role ijMadau et alJ ll9991. The 
fact that the meta galactic UV field peaks at 2 < z < 4 
ijlnoue et alJl2006l and ref. therein) points to a peak in 
the galac tic star formation (SF) rate at such redshifts 
(see also iPanter et aLll2006D . A major uncertainty in 
modeling the effect of this SF on the thermal state of the 
intergalactic medium (IGM) is the value of the escape 
fraction / esc of ionizing UV photons, which is a function 
of redshift, and galaxy type and mass. 

Best observational estimates of / csc come from the 
detailed multiwaveba nd studies of the local Universe. 
iLeitherer et alJ l)1995|) find that less than 3% of ioniz- 
ing photons from low-redshift starburst galaxies escape 
into the IGM. More recently, f or several local starburst 
galaxies. iHeckman et alJ l)2001l) estimate / esc ^ 6%, not- 
ing that inclus i on of dust will further reduce / csc , and 
iBergvall et all pOOfih find / csc ~ 4 - 10% for another 
local starburst. 

The situation is markedly different at z > 3, where 
the Lyman continuum (LyC) d etection points at much 
higher / csc ijSteidel et alJl200l|) . By comparing the di- 
rect observations of th e LyC from galaxi es to the mean 
cosmic UVB intensity, llnoue et"aTI J2006) conclude that 
the average / esc increases from 1 — 2% at z = 2 to ^ 10% 
at z > 3.6. 

Most theoretical estimates of / osc to date have been 
based on models with a smooth distribution of gas 



Electronic address: razoumov@ap.smu.ca 
Electronic address: jslarsen(9dark-cosmology.dk 

1 Institute for Computational Astrophysics, Dept. of Astron- 
omy & Physics, Saint Mary's University, Halifax, NS, B3H 3C3, 
Canada 

2 Dark Cosmology Centre, Niels Bohr Institute, University of 
Copenhagen, Juliane Maries Vej 30, DK-2100 Copenhagen, Den- 
mark 



(jR.icotti & Shull l2000j) or semi-analytical models with 
expanding shells and superbubbles in a disk galaxy 
Ipove et all 120001 iFuiita et alJ 12003). In this letter we 
present ab initio calculations of the energy-dependent 
/ csc of ionizing radiation from SF regions in proto- 
galaxies found in numerical simulations at five different 
redshifts from z = 2.39 to 3.8. We define f C sc{v : r) sim- 
ply as a fraction of photons of energy hv which reach 
a shell of radius r around a given source; by definition 
/esc( z/ , 0) = 1. We postprocess high-resolution simulation 
datasets with radiative transfer (RT) on nested grids, 
with the maximum grid resolution of 15 pc. Each galaxy 
contains from several hundred to several thousand dis- 
tinct stellar sources representing individual, young SF 
regions. Our interstellar gas distribution is extremely 
clumpy due to an array of physical processes including 
feedback from SF. We also estimate the effect of the dy- 
namical expansion of HII regions on the escape of ionizing 
photons. 

2. THE CODE AND MODELS 

To compute the escape of ionizing radiation from SF 
regions, we use an extens ion of the Fully Threaded 
Transport Engine (FTTE, iRazoumov & Cardalll 12005ft 
to point sources. This new module was previously 
briefly introduced in the comparison study bv llliev et alJ 
(2006) . It exten d s the adaptive ray-splitting scheme of 
lAbel fc Wandeld (J2002) to a model with variable grid 
resolution. Sources of radiation can be hosted by cells of 
any level of refinement, although usually in cosmological 
applications sources are sitting on the deepest level of re- 
finement (Fig. ^) . Around each point source we build a 
system of 12 x 4" _1 radial rays which split either when we 
move further away from the source, or when we enter a re- 
fined cell, and n — 1, 2, ... is the angular resolution level. 
Once a radial ray is refined it stays refined, even if we 
leave the high spatial resolution patch, until even further 
angular refinement is necessary. All ray segments are 
stored as elements of their host cells, and for actual trans- 
port we just follow these interconnected data structures 
accumulating photo-reaction number and energy rates in 
each cell. With multiple sources, we add the rates from 
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TABLE 1 

Number of distinct star-particles of age<34 Myr. 



Fig. 1. — Radial ray geometry in the FTTE point source algo- 
rithm. Color represents the angular resolution level. To avoid a 
non-uniform coverage of a cell such as that seen at the top of this 
box, we normally avoid cell interfaces with a jump in refinement 
by more than one level. 



all sources before updating the time- dependent r a te net - 
work using the chemistry solver from I Abel et al.l {1997), 
and the global timestep is chosen such that each species 
abundance does not change by more than 30%. 

With a small (N SIC < 10 3 ) number of star particles 
we can compute a full ray geometry separately for each 
source. For a large number of sources (10 3 — 10 6 ) the 
algorithm allows construction of trees of sources in which 
individual tree nodes are treated as a single source far 
from their origin. 

Our algorithm shares many common ideas with the re- 
cent RT module of the cos mological structure formation 
code Enzo implemented by lAbel et all l)2006ft to study 
feedback from Pop III stars and, in fact, is a basis of an 
independently developed radiation hydrodynamics mod- 
ule of Enzo which we will introduce elsewhere. 

2.1. Galaxy formation models 

We use results of high-resolution galaxy formation sim- 
ulations in a standard ACDM cosmology done with a 
significant ly improved version of the T reeSPH code de- 
scribed bvlSomm er-Lars en et all l)2003fl - some detail is 
given in iSommer-Larsenl l)2006T) . The simulations invoke 
the formation of discrete star "particles" , which repre- 
sent a population of stars born at the same time in ac- 
cordance with a given initial mass function (IMF). For 
the purpose of this work we focus on two (proto-) galax- 
ies, K15 and K33, which at z=0 become typical disk 
galaxies of V c = 245 and 180 kms -1 , respectively. The 
RT calculation are performed at five different redshifts, 
z = 3.8, 3.6, 3.4, 2.95, and 2.39, for "normal" resolution 
simulations, and at the first three for "high" resolution 
ones (8 times higher mass resolution, twice better force 
resolution). For galaxy K15 we also vary the strength 
of supernova feedback (Table 1). Total particle numbers 
range from - 1.6 x 10 5 to - 2.2 x 10 6 ; star (and SPH) 
particle masses are 1.1 x 10 6 and 1.4 x 10 5 M , for normal 
and high resolution simulations, respectively. 

Radiative transfer is computed on top of a nested data 
structure containing 3D distributions of physical vari- 
ables. To create such a structure from the SPH simu- 
lation datasets, for each galaxy we cut out a (250 kpc) 3 
box centered on the galaxy and projected this box onto 
a 128 3 uniform grid. We then subdivided every base grid 
cell which contains more than N max = 10 SPH particles, 
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and continued this process of subdivision recursively so 
that no cell contains more than -/V max gas particles. 

2.2. Stellar population synthesis 

In the (normal resolution) galaxy K15-1-8 at z = 3.6 
there are 347 stellar sources younger than t up — 34 Myrs 
(Fig. The birth times of particles are distributed 
almost uniformly over the span of 34 Myrs, correspond- 
ing to a nearly constant star formation rate (SFR) of 
11.6 M Q / yr. To compute the stellar UV luminosity func- 
tion for this and other galaxies, we use a population 
synthesis package Starburst 1999 l)Leitherer et al.lll999() 
with continuous SF at a constant rate 11.6M0/yr dis- 
tributed among 347 stars. For other galaxies and at other 
redshifts the SFR per particle is adjusted to account for 
the actual number of stars produced in the past 34 Myrs, 
and the specific (per unit spectrum) luminosity is dis- 
tributed uniformly among all stars in the volume. 

We assume a triple-interval Kroupa IMF with indices 
(0.5,1.2,1.7) above the mass boundaries (0.1, 0.5, 1) M Q , 
respectively, a supernova cut-off at 8 M , a black hole 
cut-off 120 M©, and solar metallicity. The resulting spec- 
trum is shown in the lower right panel of Fig. [21 

2.3. Heating and HII region expansion 

The TreeSPH code assumes the UVB of Haardt & 
Madau (1996) with self-shielding in regions in which 
the mean free path of Lyman-limit photons is less than 
1 kpc. In addition to radiative heating and cooling and 
the hydrodynamical PdV term, some heating originates 
from shock dissipation which we include in the follow- 
ing way. After we project the 3D hydro fields onto the 
grids, in each cell we apply the uniform UVB with the 
self-shielding correction and compute the amount of ad- 
ditional heating needed to keep the temperature of that 
cell constant (thermal equilibrium without sources). We 
then switch on the stellar sources retaining the extra 
heating term, and evolve simultaneously the equations 
of point source transfer and non-equilibrium chemistry 
to 10 Myrs. With continuous SF and constant UV lumi- 
nosities, all /esc reach a plateau after the first few Myrs. 
In this paper all results are presented at t = 10 Myrs at 
which point we assume convergence. 

Direct momentum transfer from UV p hotons to the 
interstellar / intergalactic gas is neglig ible ijSnitzerl 119781 
iWhalen Sz Normarill2006|) . However, heating and ioniza- 
tion by stellar photons creates pressure gradients that 
drive expansion of the HII regions into the surrounding 
gas. T he physics of expansion has been studied in deta il 
in ID ijKitavama et alJl2004i IWhalen fc Normanll20Q6|h 
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Fig. 2. — The distribution of all stellar particles (SF regions) 
younger than t up = 34Myrs in the volumes 250kpc(top left), 
60kpc(top right) and 8 kpc (lower left) on a side centered on the 
galaxy K15-1-8 at z = 3.6. The color of each particle represents 
the level of the host cell: purple circles (level 5), light blue circles 
(level 4), green circles (level 3), red circles (level 2), black circles 
(level 1), and black dots (base grid level cells). The lower right 
panel shows the input spectrum of a star particle. 



2D jShamro et al.ll20f)l . and recently 3D <|Abel et alJ 
12006(1 radiation hydrodynamics simulations. Without 
coupling RT to the hydro code, we cannot model this 
expansion self-consistently. However, we can mimic its 
effect by lowering the density in the giant HII region 
around the starburst region which will ease the escape 
of UV photons. W e use the results from section 3.1 of 
ILarsen et al.l (|200iy) which lists the " initial Stromgren ra- 
dius" Ro (1 Myr after the burst) and the "final Stromgren 
sphere" R e (4 Myr after the burst; at this time the UV 
luminosity declines rapidly) as a function of the local 
density and the ionizing source luminosity. 

For each source we use the initial gas density riff of 
the its host cell to find Rq and R e and then modify the 
density of each cell inside R e by nl rc /n^ c , provided that 
the initial density of that cell uq < rtQ rc . For uq > n.Q rc the 
expanding shell hits a denser region, and without a full 
hydrodynamical calculation one cannot predict whether 
this will lead to expansion or compression of the denser 
region, which explains our condition. Therefore, for a 
given source our density correction is anisotropic which 
seems entirely realistic for an HII front expansion in an 
inhomogeneous medium. 

For such large galaxies, we do not expect our results 
to change significantly if a fully self-consistent hydro/RT 
approach is used, as the energetic and momentum effects 
of ionization fronts are negligible compared to the effects 
of SNII explosions, shocks etc. which are already in- 
cluded into our models, resulting in a very clumpy inter- 
stellar medium (ISM). However, the radiative feedback 
would also smooth some of the small-scale structure by 
raising the Jeans mass which is especially important in 
lower mass galaxies, an effect which we plan to include 
in our future simulations. 

3. RESULTS AND DISCUSSION 
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Fig. 3. — Spectral dependence of / Gsc at r = 100 kpc for all 
five input galaxy models at (where transfer has been computed) 
z = 3.8 (red), z = 3.6 (magenta), z = 3.4 (blue), z = 2.95 (cyan), 
and 2 = 2.39 (green). For each model the redshift evolution is 
from left to right, and each curve goes from 13.6 eV to 135 eV. All 
solid lines are models without the HII region expansion, with the 
SFR distributed uniformly among all stars younger than 34 Myrs. 
For K15-06-8 the dashed lines show models with the HII region 
expansion. For K15-1-8 we also computed the models with the 
same SFR distributed among stars younger 3.4 Myrs (dotted) and 
10 Myrs (dash-dotted lines). 



Fig.|3]shows spectral dependence of / csc at r = 100 kpc 
in each model, averaged over all sources. All / esc have 
been corrected for the boundary of the volume. For in- 
dividual sources the escape of photoionizing radiation is 
similar to a phase transition: it is very sensitive to con- 
ditions at the source and therefore can vary by a sizable 
factor for the same galaxy over the course of its evolu- 
tion. We find that while there are variations among the 
models, the escape fractions tend to decrease with lower 
redshift, reflecting the fact that with time more gas cools 
into the higher density clouds, as well the declining SFR 
at z < 2.95. With normal feedback in both galaxies / esc 
at the Lyman edge drops from 8 — 10% at z = 3.8 to 
~ 1 - 2% at z = 2.39. 

Not surprisingly, in all our models / esc of photons ca- 
pable of doubly ionizing helium is very low, due to the 
softness of stellar radiation. The escape fraction of pho- 
tons capable of single He ionization is comparable at all 
redshifts to that of H-ionizing photons. 

Increasing SN feedback in K15 evacuates more gas from 
the vicinity of the SF regions leading to a threefold rise 
in /csc at z = 3.6. Lower local gas densities translate 
into a 40% reduction in the SFR (Table 1). At z = 2.95, 
a denser environment near the sources produces a much 
smaller difference in the SFR and / osc , when we increase 
the strength of feedback, and, at z = 2.39, / esc actually 
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Fig. 4. — Lyman-limit escape fraction distribution for all sources 
in galaxy K15-1-8 at z = 3.6 (without HII region expansion), at 
four different radii, from 100 pc to 100 kpc. 



drops by 25% as a larger fraction of gas swept-up by the 
shocks stays in the vicinity of the SF regions. 

Next we compare our normal and high resolution re- 
sults, at z — 3.8 — 3.4, and find a good agreement between 
the two. We see it as an additional test of our algorithm, 
since other properties of the two sets of simulations (SFR, 
etc.) also match. 

Since our choice of distributing SF over the t up = 
34 Myrs youngest stars is somewhat arbitrary, in galaxy 
K15-1-8 we also experimented with the same amount of 
SF distributed over all stars younger than t up — 10 Myrs 
and 3.4 Myrs. Ideally, these young stars would pro- 
duce most of the ionizing UV photons. Since we do 
not calculate transfer and hydro simultaneously, we feel 
that applying continuous SF with the constant rate is a 
better computational approach than assuming instanta- 
neous SF which will have no lasting effect on the thermal 
state of the IGM after few Myrs. With continuous SF we 
expect our results to converge as we increase the num- 
ber of star particles. The results at t up — 10 Myrs and 
34 Myrs are in fact very close to each other at all three 
redshifts (Fig.[3J), whereas the 3.4 Myrs results are some- 
what off, likely related to the small number of sources 
used to represent the SF in this latter case. 

The dashed lines in K15-06-8 in Fig.|31show the energy- 
dependent / csc at r = 100 kpc corrected for the HII re- 
gion expansion. At most redshifts expansion raises / esc 
by 10 — 30%, although at z = 2.95 / esc jumps by more 
than a factor of two. Further analysis shows that most 
of this change can be attributed to removal of the gas 
from the highest refinement level cells which are hosting 
the sources, not the surrounding cells. 

In Fig. 0] we show how / osc changes with distance from 
sources in galaxy K15-1-8 at z = 3.6. A large fraction 
of photons reaches the radius r = 100 pc. As we do not 



have any reliable information on gas clumping on sub- 
resolution scales at these redshifts, the effective radius 
of absorption is not firmly established at this point. In 
our current models most absorption occurs within few 
hundreds pc of each star particle, well inside the virial 
radius of ^45 kpc. At z = 3.6, the proto-galaxy consists 
of several star forming regions as well as a significant 
amount of the intergalactic HI, all of which account for 
some fractional absorption beyond 1 kpc. 

In this study we neglec ted the effect of dust extinc- 
tion. iBenson et al.l l)2006f) point out that in the case of 
a smooth ISM attenuation by dust can reduce / esc by a 
factor of few. For a clumpy ISM at high redshifts the 
magnitude of dust attenuation should be much lower, 
particularly that most ionizing photons are produced in 
star bursts, as opposed to a quiescent disk. 

In conclusion, our models with exact ionizing contin- 
uum RT around a large number of SF regions in young 
protogalaxies show a monotonic decline in / csc of ionizing 
photons from higher to lower redshifts. With the normal 
feedback strength / esc at the Lyman edge drops from 
- 6 - 10% at z = 3.6 to ~ 1 - 2% at z = 2.39. At higher 
redshifts SF on average occurs at slightly lower densities 
resulting in easier escape of UV photons into the IGM. 
T his result agrees w ell with the observational findings 
of llnoue et al.l l|2006|) on redshift evolution of / osc . Note 
that although the two galaxies in our study have a mass 
ratio of ~ 2.5, their / csc are very similar. One could sug- 
gest that perhaps, once star formation begins, it is the 
physical conditions in the clumpy gas clouds that deter- 
mine the escape of UV photons, rather than the overall 
properties of their host galaxies. However, one would 
need a much larger statistical sample to test this spec- 
ulation. In our next paper we will further examine the 
role of the fine structure of the clumpy ISM in the es- 
cape of ionizing photons, as well as the importance of 
hydrodynamical effects. 
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